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f -j-> ! Abstract 

c/3 ■ The propagation of surface (Rayleigh) waves over a rotating or- 



thorhombic crystal is studied. The crystal possesses three crystallo- 
graphic axes, normal to the symmetry planes: the half-space is cut 
along a plane normal to one of these axes, the wave travels in the di- 
rection of another, and the rotation occurs at a uniform rate about any 
of the three axes. The secular equation for the surface wave speed is 
found explicitly; in contrast to the non-rotating case, it is dispersive 
(frequency-dependent). Both Coriolis and centrifugal accelerations 
appear in the equations of motion: none can be neglected in favor of 
J^ the other, even at small rotation rates. 
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1 Introduction 

The main applications of Surface Acoustic Wave (SAW) devices are in the 
field of wireless communication, where SAW components are used to gener- 
ate and collect high frequency signals. In mobile phones, global positioning 
systems, or color television sets, these components are of course motionless. 
There exist however situations where SAWs propagate over a rotating sur- 
face. For instance, SAW sensors were used to monitor the pressure of a tire 
in a moving automobile (Pohl et al. 1997); an ultrasonic SAW-based appa- 
ratus was designed to detect flaws on the surface of a rotating member such 
as a roller bearing (Kawasaki et al. 1991); robust and shock-resistant SAW 
resonator gyroscopes have been miniaturized to a 1 cm x 1 cm surface (Jose 
et al. 2000); etc. For the theoretical aspect, Schoenberg & Censor (1973) 
seem to have been the first authors to describe Rayleigh waves in a rotating, 
homogeneous, isotropic, linear, elastic medium, for any orientation of the 
rotation axis with respect to the free surface, and with both the Coriolis and 
the centrifugal accelerations taken into account. They pointed out that, in 
contrast to the non-rotating case, bulk and surface waves are dispersive, and 
that the acoustical tensor is Hermitian instead of real symmetric; they did 
not however derive the secular equation for surface waves. Clarke & Burdess 
(1994) obtained this equation for an isotropic half-space rotating about an 
axis orthogonal to the direction of propagation and to the normal to the 
half-space. Others (Lao 1980; Wauer 1999; Grigor'evskii et al. 2000; etc.) 
considered similar problems but neglected the centrifugal acceleration. For 
anisotropic crystals, Fang et al. (2000) and Zhou & Jiang (2001) considered 
crystals of tetragonal symmetry but did not derive the secular equation ex- 
plicitly. As pointed out in the latter article, the influence of the centrifugal 
acceleration upon the SAW speed is of the same order as that of the Cori- 
olis acceleration and should not be overlooked, even at small rotation rates. 
In this paper, an analytical investigation of the propagation of SAWs in a 
rotating orthorhombic crystal is presented. 

In §2, the crystal considered is at first of general (triclinic) anisotropy 
and rotates at a uniform rate about any axis. The equations of motion 
and the boundary conditions are derived. Then the analysis is specialized 
to crystals with rhombic and higher symmetries. The material is cut along 
the plane of symmetry X2 = and the wave travels in the Ox\ direction, 
where the Ox^ axes are crystallographic axes. The semi-infinite body rotates 
along Oxi, 0x2, or Ox^. In the latter case, anti-plane motions and stresses 
decouple from their in-plane counterparts; in the other two cases, they do 
not. The choice of orthorhombic symmetry is of great importance because, as 
pointed out by Royer & Dieulesaint (1984), it covers 16 types of symmetries, 



including tetragonal, hexagonal, cubic, and of course isotropic; the choice of a 
rotation axis aligned with the normal to the free plane or to crystallographic 
axes in this plane is justified from an experimental point of view. In §3, the 
secular (dispersion) equations are derived explicitly for each rotation, using 
the method of the polarization vector. This method was used by Currie 
(1979) and by Taziev (1989) in order to derive explicit secular equations for 
SAWs in non-rotating crystals. 

Recently, Ting (2003a6) improved upon this method and placed it within 
the framework of the Stroh (1958, 1962) formalism. Here, the fundamental 
equations of the method are obtained in a manner which does not rely upon 
this formalism. The secular equation is obtained as a quadratic relationship 
between three determinants (of 2x2 matrices for rotation about Ox 3 , of 
3x3 matrices for rotations about Ox\ or Oxz), whose elements are given 
explicitly in terms of the elastic parameters, the wave speed, and 5 - - the 
ratio of the rotation rate by the wave frequency. Several secular equations 
arise in the resolution of the problem but only one is kept; the others are 
dismissed on the argument that they do not reduce to the known single cubic 
secular equation when the rotation rate vanishes. Numerically, the Rayleigh 
wave speed decreases monotonically as 5 increases. This influence is shown 
graphically for a- Iodic Acid (HIO3, rhombic). For this graph, 5 is chosen to 
vary from to 10. Although it is not realistic to imagine a crystal rotating 
at a frequency which would be greater, or even close to, the frequency of 
an ultrasonic surface wave, this wide range is chosen to show numerically 
that the wave speed tends to zero with increasing S without, seemingly, ever 
reaching it. For Silica (SiC-2, isotropic), the calculations are conducted for 
a smaller range (0.0 < S < 0.1) with and without the contribution of the 
centrifugal acceleration in order to show its importance. 

Some theoretical issues are left open such as, the questions of existence 
and uniqueness of a surface wave in a rotating crystal, or the question of the 
general behaviour of the wave speed as a function of 5. 

2 Equations of motion 

2.1 General anisotropy and arbitrary rotation axis 

Consider a semi-infinite body made of a linearly elastic, anisotropic crystal. 
For two-dimensional deformations, the stress-strain relations are (Ting 1996, 
p. 38), 

<x° = C°e°, or e° = s'<x°; CV = I, (1) 



where 
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Here the CV,- are derived, using the Voigt (1910, p. 560) contracted notation, 
from the fourth-order elastic stiffness tensor and the reduced elastic compli- 
ances s'^ are defined by ^3. 

In the Cartesian orthogonal coordinate system (Ox\x 2 x 3 ) = (Oijk) where 
x 2 > is the region occupied by the half-space, a surface (Rayleigh) wave 
traveling at speed v and wave number k in the xi-direction with attenuation 
in the ^-direction is described by the following mechanical displacement 
field, 

u(x u x 2 ,x 3 ,t) = U(kx 2 )e ik(xi - Vt \ U(oo) = 0. (4) 

From ([Tj)i it follows that the corresponding stress components are of similar 
form. In particular, the tractions acting upon the planes x 2 = const can be 
written as 

a l2 ( Xl ,x 2 ,x 3 ,t) = ikt i {kx 2 )e lk ^- vt \ (5) 

Now consider that the crystal is rotating at a uniform rate Q about the di- 
rection of a unit vector w. Then the incremental (time-dependent) equations 
of motion in the rotating frame are (Schoenberg & Censor 1973) 



div cr = pii + 2pfiw x u + pfi 2 w x (w x u) 



(6) 



where p is the mass density and the dot denotes differentiation with respect 
to time. Note that on the right hand-side of the equations, both the Coriolis 
(second term) and the centripetal (third term) accelerations appear. Using 
(PQ) 5 (SD , and ([S]), these equations can be formulated as a linear homogeneous 
system of first-order differential equations for U = [Ui,U 2 ,U 3 ] T and t = 
[ti,t 2 ,t 3 ] : 
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where X ■ 
Here, N 1; 



pv 2 and the prime denotes differentiation with respect to kx 2 . 



N 5 



N 2 T , and N 3 



N 3 T are the usual submatrices of the 



fundamental matrix N (Ingebrigsten & Tonning 1969). Barnett & Chadwick 
(1990) and Ting (1988) give explicit expressions of the components of the N; 
in terms of the CV, and of the s^ -, respectively. The matrix J contains all the 
information relative to the rotation. Explicitly (Schoenberg & Censor 1973), 



Jij = [(Sij - WiWj)5 2 - 2ie ijk WkS}X, 5 = Q/u, 



(8) 



where u = kv is the real frequency of the wave, 5^ is the Kronecker operator, 
and €ijk is the alternator. Note that J is Hermitian. 

Finally, the following boundary conditions apply: vanishing of the trac- 
tions on the free plane surface X2 = and of the wave as x<i — > oo, that 
is 

t(0) = 0, U(oo)=0, t(oo) = 0. (9) 

2.2 Orthorhombic crystals 

Henceforward, we are concerned with orthorhombic crystals whose crystallo- 
graphic axes are aligned with the Oxi. In that case, we have 



Ni 






-1 


o" 


; '12 








4i 
















N 5 



°66 










S 22 



b L2 



s 






/ 

44 



(10) 



Also, the crystal is assumed to be rotating about a crystallographic axis. 
Introducing the Hermitian tensor IC 1 ) = N3 + J + XI (lower left submatrix 
in dZD); we have 
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for rotation about the Ox\ axis, OX2 axis, and Ox^ axis, respectively. 



3 Explicit dispersion equations. 

The Stroh-Barnett-Lothe formalism has beautifully addressed and resolved 
most of the problems arising in linear anisotropic elasticity (see the textbook 
by Ting (1996) for a thorough review.) Some researchers have however fol- 
lowed other routes for the study of SAWs in anisotropic crystals (see Ting 
(20036) for an exposition of several methods.) For instance, Fu & Mielke 
(2002) recently proposed a highly efficient numerical scheme to determine 
the Rayleigh wave speed, which is based on a Ricatti equation. Here the 
method of the polarization vector (Currie 1979; Taziev 1989; Ting 2003ao) is 
presented for rotating orthorhombic crystals. Only simple algebraic manip- 
ulations and first integrals (Mozhaev 1995) are used, so that the derivation 
should appeal to the reader who is not familiar with the Stroh (1958, 1962) 
formalism and with the orthogonality relations (Barnett & Lothe 1973). 

3.1 Method of the polarization vector. 

Let us examine the equations of motion (J7|) closer. They read: 

- iU' = NiU + N 2 t, -it' = K (1) U + Ni T t. (12) 

Take the scalar product of the first line by t and of the second line by U ; 
then add, together with the complex conjugate; finally, integrate between 
and oo to obtain: 

[U • K (1) U + t ■ N 2 t + t ■ NxU + t • NiU]g° = 0, (13) 

which, using the boundary conditions Q), reduces to: U(0) • K^TJ^O) = 0, 
a relation first established by Stroh (1958) in the static, non-rotating case. 

This procedure is easily generalized as follows. First compute N (the 6x6 
matrix on the right hand-side of the equations of motion ([7|)) to the power 
n, where n is any positive or negative integer, as 
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(14) 



(say). Of course, the matrices N-j 11 , N 2 , and K^ n ^ do not coincide with 
Ni™, N2™, and K n , respectively. The matrix Nij is however symmetric just 
like N 2 is, as can be proved by induction (Ting 2003a). Similarly it can be 
proved, in essentially the same manner, that IC n ) is Hermitian just like JO 1 * 1 



is. Now multiply (j7|) in turn by N l and by — iN to get 

K-i)tt/ , tm(-i)+/ r -iNxU' - iN 2 t' = Ni 2) U + N^ 2) t, 



iU = Nf XJ U' + NipV, 
it = K^U' + Ni" 1)T t', 



and 



-iK«U' - iN! T t' = K^U + N< a) t, 

(15) 
Take the scalar product of the first line in the first (second) system by t 
(t ) and of the second line by U (U); then add, together with the complex 
conjugate; finally, integrate between and oo to obtain: 



[U • K (_1) U + t • N 2 _1) t + t ■ N { ~ 1} V + t ■ N[~ 1} V} 
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r( 2 )] 



[U • K (2) U + t • N 2 2) t + t • N^U + t • Ni 2, U]S° = 0, (16) 

respectively. Using the boundary conditions flOJ), we have: U(0) ■K^ 1 ^U(0) = 
and U(0) • K^ 2 ^U(0) = 0. These steps may be repeated ad infinitum for 
any positive or negative power n of N to give 



U(0) • K (n) U(0) = 0. 



;i?) 



However, because of the Cayley-Hamilton theorem, only five N n are linearly 
independent and consequently, (IT7|) yields at most five equations. In a non- 
rotating frame, Currie (1979) used the equations written at n = 1,2,3 for 
SAWs in the plane of symmetry of a monoclinic crystal and Taziev (1989) 
used the equations written at n = 1, 2, 3, 4, 5 for triclinic crystals. Recently, 
Ting (2003a) pointed out that the choices n = —2, —1, 1, 2, 3 lead to simpler 
expressions. 



3.2 Rotation about Ox\ (propagation direction). 

When the orthorhombic crystal rotates about the propagation direction Ox\. 
the matrices K^ 2 ^ and K^ 2 ^ have the following structure, 
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where K^ and K^ are real. Explicitly, 
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and 
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Equations (|T7|) written at n = —2, 2, for U(0) = [1, a± + ia 2 , Pi + i/32] T give: 
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Hence, either (a) a{~ 2) a{ 2) - k[fK^ 2) = or (b) ai = fo = 0- The first 
possibility is not a valid secular equation, because it is not consistent with 
the non-rotating case; in other words, (a) at 5 = does not reduce to the 
cubic Rayleigh function for non-rotating orthorhombic crystals, 

1 - (2s' u + s' m )X + s' n {s' u - s' 22 + 2s' m )X 2 + s' n (s' u s 22 - s'? 2 - s' n s' m )X 3 = 0. 

(22) 
Thus (b) applies and 

U(0) = [l,ia 2 ,/3 1 ] T . (23) 

The directions of the vectors i + /3]k and a 2 j are conjugate directions with 
respect to the polarization ellipse of the wave at the free surface (Boulanger 
& Hayes 1993). In fact, because they are orthogonal, they are along the 
principal axes of the ellipse. Now we compute a 2 and 0±. 

The matrices K^ _1 \ IC 1 ), and K^ 3 ^ have the following structure, 
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where Kj^ , fQ 2 , fQg , and -fQg are reoi. Explicitly, 
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(where D\ is defined in (I201U ); the components of YS^ are given in (ITTjb ; 
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Substitution of (1251) into equations (fTTl) written at n = —1,1,3 gives: 
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These three equations are rewritten as F^gu = —hi with 
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Let A = det F, and A& be the determinant of the matrix obtained from 
F by replacing the fc-th column with h. Then the solution to Fg = h is 
g k = -A fe /A. But g\ = 4gig 2 , that is 

A 2 - 4A!A 2 = 0, (29) 

the explicit dispersion equation for Rayleigh waves in rhombic crystals rotat- 
ing about Ox\. It is a polynomial of degree 12 in X = pv 2 and of degree 10 
in S 2 . At S = 0, A23 = 0, so that Ai = A 2 = 0, while A3 factorizes into the 
product of a quadratic in X and the cubic (122]) . 



3.3 Rotation about Ox^ (normal to the free plane). 

When the orthorhombic crystal rotates about Ox 2 , the normal to the free 
plane, the matrices K4 2 ) and IC 2 ) have the following structure, 



K(n) 





k {n) 
-"12 





-"12 









\K. 



(n) 



-iA. 



23 



23 







n = -2,2, 



(30) 



where AJ2 an d -^23 are rea ^- Explicitly, 
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0, or (b) «! = a 2 = 0, or (c) 
Oi\ = Pi = 0. The first possibility is not a valid secular equation, because it is 
not consistent with the non-rotating case (that is (a) does not reduce to the 
cubic ( 122]) at 5 = 0.) It can be checked that in Case (b), the equations ( TTT)) 
written at n = —1, 1, 3, lead to a secular equation which vanishes identically 
at 5 = 0, so that (b) is also dismissed. Thus (c) remains, and 
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The directions of the vectors a 2 j + /3 2 k and i are along the principal axes of 
the polarization ellipse at the free surface. 

The matrices K> -1 ), K^ 1 ), and JC 3 -* have the following structure, 
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where K\ x , K 22 , K^ 3 , and K\ 3 are real. Explicitly, 
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Substitution of (1511) into equations (TT71) written at n 
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These three equations are rewritten as F^g^ = —hi with 
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Let A = det F, and A& be the determinant of the matrix obtained from 
F by replacing the fc-th column with h. Then the solution to Fg = h is 
-Afc/A. But g 2 = 4g 1 , that is 



9 k 



A 2 + 4AxA = 0, 



(40) 



the explicit dispersion equation for Rayleigh waves in rhombic crystals rotat- 
ing about Ox 2 . It is a polynomial of degree 12 in X = pv 2 and of degree 10 
in S 2 . At 5 = 0, K]^ = 0, so that A = A x = 0, while A 3 factorizes into the 
product of a quadratic in X and the cubic (T22|) . 

Figure 1 shows the influence of 5 upon the surface wave speed for a- 
Iodic Acid (HIO3, rhombic). The elastic stiffnesses (10 10 N/m 2 ) are (Royer 
& Dieulesaint 1996): C n = 3.01, C 12 = 1.61, C 22 = 5.80, C 44 = 1.69, 
C 55 = 2.06, and C 66 = 1.58; the mass density is 4640 kg/m 3 . We see that 
the Rayleigh wave speed is a monotone decreasing function of 5. When the 
crystal rotates about the X2-axis (solid top curve), the speed is greater at a 
given 5 than when it rotates about the Xi-axis (solid bottom curve) or the 
23-axis (dotted curve). 
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Figure 1: Rayleigh wave speeds for a-Iodic Acid (rhombic) rotating about 
Ox\ (bottom solid curve), 0^2 (top solid curve), or Ox 3 (dotted curve). 

3.4 Rotation about Ox?,. 

The case where the orthorhombic crystal rotates about an axis orthogonal to 
the direction of propagation and to the normal to the free surface has been 
treated elsewhere (Destrade 2003) as a special case of a monoclinic crystal. 
It turns out that the in-plane strain decouples from anti-plane strain, and 
that the displacement is of the form, 



U(0) = [l,ia 2 ,0] r 



(41) 



where 0:2 is real. Explicitly, the wave propagates at a speed which is a root 
of the following dispersion equation: 
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with 
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F 22 = -2s' u 8X, F 23 = s' n (l + 5 2 )X, 

h = [s 22 (l + 5 2 ) - (s' lv s' 22 - 4)(1 - 5 2 ) 2 X]X, h 2 = l- s' u (l + 5 2 )X, 

(43) 
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and the mechanical displacement at the free surface is given by (|41l) , with 
«2 = \F12h2 — F 2 2h 1 ) / (Fi 2 F 23 — F 22 F 13 ). 

The dispersion equation is a polynomial of degree 6 in X = pv 2 and 
in 5 2 . In the isotropic case, it reduces to a^ 1 = (i = 0, 1, . . . , 6), where 
£ = pv 2 /C 66 , a = C m /Cn, and 

o 6 = (1 - 5 2 ) 4 [(l + «) 2 (1 + 5 2 ) 2 - 4a(l - 5 2 ) 2 ], 

a 5 = -16(1 - a)(l + 5 2 )(1 - 5 2 ) 2 [(l - a)(l + <5 4 ) + 2(1 + 3a)<J 2 ], 

a 4 = 16(1 - a) 2 (l + 5 2 ) 2 [(7 - 2a)(l + 5 4 ) + 2(1 + 2a)<J 2 ], 

a 3 = -32(1 - «) 2 (1 + <5 2 )[(13 - 9a)(l + <5 4 ) + 2(7 - lla)<5 2 ], (44) 

a 2 = 64(1 - a) 2 [(13 - 16a + 4a 2 )(l + <5 2 ) 2 - A5\ 

ai = -256(3 - 2a)(l - a) 3 (l + 5 2 ), 

a = 256(l-a) 4 . 

This equation corresponds to the rationalized form of the dispersion equation 
obtained by Clarke & Burdess (1994). In the non- rotating case, 5 = and 
the sextic secular equation reduces to the (squared) cubic of Rayleigh (1885): 
£3 _ g£2 _|_ ^4 _ 16a)£ — 16(1 — a) = 0. Finally in the incompressible isotropic 
case, a = and the sextic equation factorizes into the product of the two 
following cubics, 

(1 + 5 2 )(1 ± 5) 4 £ 3 - 8(1 + <5 2 )(1 ± 5) 2 ¥ + 8(3 ±25 + 35 2 )£ - 16 = 0. (45) 

Figure 2 shows the influence of 5 upon the wave speed (thick curve) for 
Silica (SiC>2, isotropic) rotating about the normal to the directions of prop- 
agation and of attenuation. The elastic stiffnesses (10 10 N/m 2 ) are (Royer 
& Dieulesaint 1996): C\\ = 7.85 and C\ 2 = 1.61; the mass density is 2203 
kg/m 3 . The thin curve corresponds to the dispersion curve obtained when 
the centrifugal acceleration is omitted in the equations of motion (JO}. We 
see that even for small 5, this term plays an important role: for instance, the 
speed is reduced from 3409 m/s in the non-rotating case to 3350 m/s (hori- 
zontal dotted line) for 5 = 0.08852 with the full equations of motion, and for 
5 = 0.04878 without the centripetal acceleration, a relative difference of 45% 
for the expected value of 5\ Similarly, Grigor'evskii et al. (2000) argued that 
one of the bulk waves disappeared in a rotating isotropic media at: 5 = |; 
in fact, as shown by Schoenberg & Censor (1973), this phenomenon occurs 
at 'resonance': 5 = Q/u = 1. 
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Figure 2: Rayleigh wave speeds for Silica (isotropic) rotating about Ox% 
calculated with (thick curve) and without (thin curve) the contribution of 
the centrifugal force. 
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3.5 Remark: rotation about any direction in x% = 
plane. 

Note that the method of the polarization vector may also be carried out fully 
in the case where the crystal rotates about any axis in the {x\X 2 ) plane. Then 
the secular equation is obtained as a nonlinear combination of determinants 
of 5 x 5 matrices, as follows. 

When the unit vector w along the axis of rotation is such that w% = 
and W\W 2 7^ 0, then the matrices IO n ) have the following structure, 



K(n) 



21 11 -^12 lrv 13 

kg kg wg 

-iKg -iKg kg 



(46) 



where kg, kg , kg, kg, Kg, and Kg are real. 

Substitution of U(0) = [1, a.\ + ia 2 , /?i + i/?2] T into equations ffTTj) written 
at five different n gives: 

7->( n ) 1 r>( n ) 1 r>( n ) 1 r^( n ) T^( n ) r>( n ) 

#22 0i + K 33 'g 2 + K{ 2 >g 3 + K{ 3 > g 4 + Ky g 5 = -K\ x ' , 

g l=a \ + a\, g 2 = ^ + /3 2 , g 3 = 2a u g 4 = -2/3%, g 5 = 2(a 2 /3 1 - «i/3 2 )- 

(47) 

The nonhomogeneous linear system of five equations is solved for the five 
unknowns g^. But these quantities are related through 

g\gi + g\gi + g\ - gig^gs - % 2 ^4 = 0, (48) 

the dispersion equation for Rayleigh waves in rhombic crystals rotating about 
any axis in the (x\x 2 ) plane. 

Explicitly, the components of the IC n ) matrices are however too long 
to reproduce here. Also, the corresponding set-ups are unlikely to be of 
experimental relevance. 
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